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Abstract. - We consider the prototypical "piston pump" operating on a ring, where a circulating 
current is induced by means of an AC driving. This can be regarded as a generalized Fermi-Ulam 
model, incorporating a finite-height moving wall (piston) and non trivial topology (ring). The 
amount of particles transported per cycle is determined by a layered structure of phase-space. 
Each layer is characterized by a different drift velocity. We discuss the differences compared with 
the adiabatic and Boltzmann pictures, and highlight the significance of the "diabatic" contribution 
that might lead to a counter-stirring effect. 



Stirring is the operation of inducing a DC circulating 
current by means of AC driving. This is naturally achieved 
by integrating a pump [Ifj3] in a closed circuit [4j[5]- 
It can also be regarded as a variation of a Hamiltonian 
ratchet plm where transport is induced in a periodic ar- 
ray. Pumping and stirring have largely been considered 
in the regime of slow (adiabatic) driving, where it can be 
related to the Berry phase that is associated with the driv- 
ing cycle. This adiabatic approach is based on a simple 
picture of probability flow. 

Challenging this oversimplified view, we argue that 
there are typical circumstances where the analysis should 
go beyond the adiabatic picture, even for very slow driv- 
ing. We here present a detailed account of determinis- 
tic stirring that naturally extends into the non-adiabatic 
regime, complementary to related studies of stochastic 
stirring |8 , Brownian ratchets 9 , Brownian motors 
10 , stochastic 11 and chaotic [12] pumps. 

that for a prototype system, the 



stochastic 11 

We shall show tnat tor a 
oscillating-piston model, even if the driving is very slow, 
the dynamics is actually complex, due to a non-trivial 
structure of phase-space, leading to drastic consequences 
for the transport. 

Outline. — After introducing the model, we describe 
the expectations that are based on a stochastic Boltzmann 
picture, and on a deterministic adiabatic picture. These 
suggest two different parametric results for the amount Q 
of pumped particles. Then we present a proper analy- 
sis of the mixed phase-space dynamics, and highlight the 
limitations of the traditional reasoning. 
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Fig. 1: (Color on-line) An oscillating piston pump is integrated 
into a ring. The pumping cycle is illustrated on the right. The 
scattering within the ring is modeled as a potential barrier. 
We are considering a closed geometry that can be regarded 
a generalization of the Fermi-Ulam model. If we considered 
an open geometry, the pumping device would be connected 
between two unbiased reservoirs (not presented). 



The model. — We consider the prototype system il- 
lustrated in Fig.[T} A particle with mass m moves in a 
ring of length L. There is a fixed barrier Vb at x=0 and a 
moving wall (piston) of height V p at x=X p . In the oscillat- 
ing piston paradigm |3}[l3] , the two control parameters of 
the piston are cycled periodically through a closed loop in 
parameter space. The pumping cycle (X p (p), V p (tp)), with 
(f = Qt, consists of translating the piston some distance to 
the right, shrinking its "height", pulling it back to the left, 
and restoring its original "height" . In the sequel, we shall 
assume a harmonic driving with phase shift tt/2 between 
the two parameters, 



X p ((p) = X - SX P cos((/?), 
V p ((p) = V p + SV P sin(<p). 



(1) 
(2) 



The velocity v of the particle in dimensionless units is 
defined as u = v/(SX p Q). With the same convention the 
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piston velocity X p is expressed as u x ((p) = sm(cp). The 
transmissions of the barrier and the piston are given by 
boolean expressions (true=l, false=0): 



9b(u) 



\u\ > u h 



\u-u x ((p)\ > U p ((f) 



(3) 
(4) 



where u^^ = [2T4,p/ m ] /(SX P Q). In addition to the 
three dimensionless parameters (ui),u p ^Su p ) that describe 
the barrier and the piston, we specify the geometry 
of the system defining 2n£+ = Xo/5X p and 2tt£- = 
(L— Xq)/SX p . A Poincare section of the dynamics is ob- 
tained by taking snapshots of ((/?, u) after each collision 
with the fixed barrier: 



27r£ n 



(5) 



-u n + 2u x (<p' n ), if g p (u n ,(p f n )=0, (6) 



if^«)=0. 



(8) 



Above £ n (£ f n ) is the scaled travel distance from the bar- 
rier to the piston (from the piston to the barrier). In the 
"static-wall approximation" it is £+ or £_ depending on 
the sign of u ni while in the simulations the exact value 
can be numerically determined. 

Objective. — The map above generalizes the Fermi- 
Ulam model (FUM) |14 : here we have finite heights 
of piston and barrier, and periodic boundary conditions. 
Furthermore, while the FUM has been conceived to study 
energy absorption due to deterministic diffusion in mo- 
mentum, here our interest is in the directed transport 
along the spatial coordinate. The amount of particles that 
are pumped per cycle is: 
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(9) 



where the current / = pNv((p) at a given moment of time 
is expressed by the spatial density of the particles pN, 
and the drift velocity v((p). Assuming that a first-order 
description of adiabatic transport applies, if the piston is 
displaced with velocity X p , one expects a well defined in- 
duced current / = Tlp^X pi where 7Z is an (^-independent 
dimensionless coefficient that primarily depends on V p . 
One expects 7Z — >• 1 for V p — » oo and < 7Z < 1 for fi- 
nite V p . In such a case the amount of particles that are 
pumped per cycle becomes a parametric integral: 



v(<p) = K((p) X p 



Q = p N j> K(<p)dX p . (10) 



Below we try to formulate an adequate picture of the dy- 
namics that interpolates between the opposite extremes of 
exclusively chaotic and purely regular motion. In partic- 
ular: (i) we clarify what is the drift velocity for a general 



non-equilibrium steady state; (ii) we discuss whether a 
first-order description of adiabatic transport applies. We 
first review the common expectations. 

The stochastic Boltzmann picture. — The piston 
stirring problem has been analyzed in the past 



13 



using 

a stochastic approach with transmissions < gb,g P < 1. 
Motivated by the prevailing literature that focuses on elec- 
tronic systems one assumes a Fermi-like energy distribu- 
tion, such that the initial occupation is 



f(x,p) 
Pn 



1 

rrm F 
irH 



for \p\ < mv F (11) 
[density of particles] (12) 



where mv F is the Fermi momentum. In a classical con- 
text, h can be regarded as a parameter that determines 
the occupation density. Of interest is f(p) = /(0,p), the 
momentum distribution at a section x = through which 
the current is measured. Its integral over dp gives the 
density of particles pN- If \X p \ <C v F , the momentum ex- 
change due to collisions with the piston affects only a nar- 
row shell of width ~ 2mX p around the Fermi energy E F . 
Accordingly one writes 



dp[/(p) - K-pMp) = pn nx p 



(13) 



From here Eq. (10) is implied, while 1Z is related to the 
reflectivity of the piston as follows [l3]: 



(1 ~ gg) 9b 
g p +g b - 2g p g b ' 



(14) 



In the absence of a barrier (g^ = 1) the result is 7Z = 1. 
If the ring is like a reservoir = 1/2) one observes that 
1Z = l—g p is the reflectivity of the piston. The latter result 
is very well known [3], conventionally derived using the 
scattering formalism [TJ[2] . 

Deterministic adiabatic picture. — In our model 
the transmission coefficients of the barrier are given by 
Eqs.(|3]|4|. In Eq. ([l4| we have to substitute the val- 
ues at the Fermi energy E F = (1/2) m^, which gives 
TZ((p) = {Oorl} depending on whether the piston is "be- 
low" or "above" the Fermi energy. In the latter case, from 
Eq.dlOk with Eq.Q, we get dQ = [mv F /7rh]dX p , which 



coincides with [3] . If the Fermi energy is "above" the pis- 
ton during the whole cycle, meaning that V p (cp) < E F for 
any (/?, we find Q = 0. However, in the latter case, there is 
a naive (wrong) picture that suggests a finite result: As- 
suming that Q is determined by the fraction of particles 
that are affected by the motion of the piston, the effective 
Fermi energy is V p (<p), and hence 



1 



2mV p ((p) dX p . 



(15) 



This "area" that is enclosed by the cycle resembles the ac- 
tion integral that is encountered in the canonical adiabatic 
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Fig. 2: The (ip, u) Poincare section for a particle in a ring, with /±=15, considering several cases: (a) ring with u p =oc and i£b=44; 
(b) here Ub=0, and u p =26, and (5t£ p =2.04; (c) here Ub=§, and u p = 380 and 5u p =2.55\ (d) the same but with Ub=378. The color 
of each trajectory reflects its average drift velocity (u). Non transporting trajectories are colored gray. The dashed red lines 
indicate the piston scattering threshold it ± (cp), while the dashed blue line gives that of the barrier Ub. The dot-dashed black 
lines are calculated adiabatic trajectories. Above the barrier the variation in \u\, due to the bounces with the piston, is sin(<^), 
with alternating sign in each bounce. Below the barrier u keeps its sign, and the adiabatic variation is [(|ifc|)/(27rZ+)] cos(y>). 
See the appendix for details regarding the adiabatic trajectories. 



picture. It has the form of Eq. ( |Tq| ) with a modified 're- 
flection' coefficient TZ((p) = [V p (ip) /E^] 1 ^ 2 . In spite of the 
wrong reasoning, Eq.(Jl5| is interesting because it can be 
justified as an approximation to what we call later "adia- 
batic contribution". 

Non-adiabatic deterministic dynamics. — The 

Poincare section for the generalized FUM Eqs.([5][8| is il- 
lustrated in Fig[2] We indicate there (by dashed red lines) 
the threshold velocity for piston reflection 

uf((p) = u x (ip)±u p (ip) (16) 

which is implied by Eq.Q. We define as the veloci- 
ties that correspond to max[w+] and min[iz~], respectively, 
and denote by the associated kinetic energies. For 

simplicity of presentation we assume that V& is smaller 
than E^~\ which is always smaller than E^\ Accord- 



ingly, the ballistic motion of clockwise moving particle is 
not affected if E > E^~\ while for anticlockwise moving 
particles the condition is E > E^\ 

The non-integrable region in phase-space occupies the 
rectangular strip [v^~\v^]. One expects adiabatic dy- 
namics if the slowness condition \u\ ^> £ is satisfied there. 
Looking at Fig^d one observes that this region consists of 
layers. In each layer the motion is chaotic, with the pos- 
sible exception of regular motion in some small islands. 
The lower layers, labeled collectively as r = 0, are non- 
transporting: either the particle is bounded in the £+ or in 
the £- segment, or else it occupies the whole ring without 
being ever able to cross one of the barriers. The motion in 
the subsequent transporting layers, labeled r = 1,2, ...r c , 
is characterized by a non-zero drift velocity 

(v) r = u r SX P (17) 
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Stirring. — Define as the minimal energy re- 

quired for transporting motion. Consider a zero temper- 
ature Fermi occupation as the preparation. If either E F 
is below or above E^ the induced current would 
be zero. In the latter case the non-regular transporting 
motion is merely a "musical chair" dynamics that takes 
place deep in the Fermi sea. The current would be non- 
zero if E^ < E F < E^\ meaning that only a fraction of 
the transporting region is occupied. The amount of parti- 
cles that are pumped during a period is Q = (2tt/Q)pn(v} , 
where the cycle- averaged drift velocity is 



(18) 



The normalized occupations satisfy ^2 f r = 1, where 
r = 0, 1, 2, r c . For a saturated occupation the f r are 
proportional to the phase-space area of the filled layers. 
In the latter case (v) is merely the average velocity within 
the occupied region. If the whole transporting region is 
saturated, one obtains 



(v) 



JO- 



,(+)! 



,(-) 



(19) 



It should be realized that the dimensionless amplitude of 
the piston velocity is unity. Accordingly the condition for 
the emergence of multiple transporting layers is Su p > 2. 
Numerical results for (v) r are presented in Fig.[3]3, and 
additional plots are available in Fig.[4j Typically (v) r is 
dominated by the the contribution that comes from the 
free ballistic stage of the motion. This contribution can 
be estimated quite easily (see below). 

The drift velocity. — At this stage we have to clar- 
ify how the drift velocities (v) r are determined by the dy- 
namics. The winding number (WN) of a trajectory that is 
generated by Eqs.([5][8|, and its duration, are respectively 

WN = ^2g h (v n )sign(v n ), Time = ^ (20) 



The drift velocity is 



WN _ Y, n 9b(v n )sign(v n ) 



Time 



g b (v n )sign(v n ) / \vn\~ 1 



(21) 



Using ergodicity the time average can be replaced by 
phase-space average. We just have to remember that cp 
is like time, which is canonically conjugate to the energy. 
Hence the integration measure in phase- space is ex vdvdcp. 
Consequently we get the obvious result 



v = 



JJ g b (v) v dvd(p = 
JJ dvdif 



dip 
2tt 



(22) 



where the phase-space integration extends over the layer 
that is filled by the trajectory, and v(ip) is defined as the 
parametric drift velocity. Results for the drift velocity are 
presented in Figl3k. 
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Fig. 3: (Upper panel) The drift velocity u{<p>) for trajec- 
tories in 4 representative layers (blue, cyan, orange, green) of 
Figj2ji. The motion consists of "free" and "adiabatic" stages. 
If the strict adiabatic picture were applicable it would be de- 
scribed by an expression 7Z((p)X p ((p), with < H(ip) < 1, as 
illustrated by the dotted curve. The dashed curve is X(ip). 
The dashed horizontal lines correspond to Note that 

the [—2, +2] range of the vertical axis is zoomed. (Lower 
panel) The drift velocity (u) as function of (\u\), for the var- 
ious layers that appear in Fig|2ji (diamonds). One should ex- 
clude the framed data which represent additional layers that 
appear in Fig[5]: (circles). The dashed black line assumes u > 
motion during the ballistic stage (see appendix), while the ver- 
tical dotted lines correspond to values of u that accommodate 
an even number of bounces. 



Estimate of the 'free' contribution. — (can be 
skipped in first reading). Looking in Fig.]^ we observe 
that in the absence of a barrier the motion in the trans- 
porting layers is not strictly adiabatic. Rather it contains 
a "diabatic transition" from adiabatic to free motion when 
g p switches from to 1. For u > this transition takes 
place at the intersection of the stochastic layer with u + ((p). 
The relative filling of the u < branch is determined by 
the number of bounces that are accommodated in the adi- 
abatic stage. It follows that the dependence of the drift 
velocity on (\u\) is modulated as seen in Figj3]o- 

Neglecting small uncertainty of order X/v, the ballistic 
motion takes place in the Acp interval where \u\ > u p ((p). 
Hence this interval is determined by the roots of the equa- 
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Fig. 4: Additional Poincare sections (top) showing mixed 
phase space. There is no barrier. The parameters are l± — 15, 
while u p = 27 (left) and u p = 29 (right) with Su p = 2.24, 2.48, 
respectively. The color code and lines are as in Fig[2]D. Corre- 
sponding plots of (u) in the bottom panels are as in Fig[3] 

tion Up + Su p sm(cp) = \u\. Assuming that the ballistic 
motion takes place in the u > region of phase- space, 
without branching to the u < region, we get the upper 
bound estimate 



<«>r = <H>('^) 



(23) 



An improved estimate should take the branching into ac- 
count. This branching depends on the number of bounces 
that are accommodated in the adiabatic stage, i.e. dur- 
ing the interval Acp c = 2ir — Acp. The variation of cp be- 
tween successive bounces and the number of accommo- 
dated bounces are respectively 



Acp 



2tt£ 

(MP 



A(f c 
A(po 



(24) 



We expect the branching to become negligible each time 
that m crosses an even integer number. This in confirmed 
by the numerical results of Figj3}}. 

Discussion of the adiabatic picture. — Canonical 
adiabatic theory [l5] ensures the invariance of the actions 
for a sufficiently slow perturbation of an otherwise inte- 
grable system. It is analogous to the conservation of the 
energy-level-index in its quantum version [l6] . It holds as 
long as the trajectory does not change its topology. If a 
trajectory does change its topology, we call this occurrence 
a diabatic transition. 

In the oscillating-piston model, each time that the mo- 
tion switches between free ballistic motion, extended adi- 
abatic motion, and bounded adiabatic motion, it is a dia- 
batic transition. The two types of adiabatic motions are 
described in the appendix. 

The small parameter of the adiabatic theory is X. In the 
absence of magnetic field the zero-order adiabatic states 
carry no current. This is true both classically and quan- 
tum mechanically: note that in the latter case the para- 
metric eigen- function \X) is real. It follows that adiabatic 



transport requires to go beyond zero-order. Linear re- 
sponse theory is based on a first- order treatment, leading 
to the Kubo formula for the transport coefficient. In open 
geometry the scattering formalism leads to the same re- 
sult. 

Using a quantum language, but referring on equal foot- 
ing to the classical picture, the evolving zero-order adia- 
batic state \X(t)) does not satisfy the continuity equation: 
at any moment (/) =0. Still we can deduce from the zero- 
order parametric solution a non-zero result for the current. 
This is done by associating a parametric velocity to each 
"piece" of the evolving probability distribution. This leads 



to Eq.(15). We note that such procedure has been used in 



|17| . We also note that such procedure becomes ambigu- 
ous in the case of multiple path geometry: to associate a 
"displacement velocity" to each piece in phase space is not 
always well defined. 

Contrasting with adiabaticity. — The observed re- 
sults imply that even for very slow driving the analysis 
should go beyond the adiabatic picture. If a first-order 
adiabatic transport picture were applicable, the drift ve- 
locity of the particle would adjust itself to the motion of 



the piston, and it would be possible to write Eq.(10) with 



an ^-independent 1Z. In practice we observe that (v) r 
has no simple linear relation to X p . The deterministic 
adiabatic result Eq. ( 15 ) would be obtained if we had a 



saturated occupation of the region bounded by Up(ip). It 
can be regarded as a rough approximation to the total 
adiabatic contribution that we discuss below. 

We already pointed out that the value of (v) r within 
a layer of phase space is a sum of adiabatic-like and free 
ballistic contributions. The adiabatic-like contribution is 
proportional to the integral over X p during the g p = 
stage. The free ballistic contribution is due to the average 
velocity in the g p = 1 stage, during which the particle cir- 
culates the ring without being back-scattered. The result 
of the free ballistic contribution depends on the branch- 
ing that has been explained previously, and accordingly, 
due to the alternating branching ratio, we observe an al- 
ternating net result for (v) r as we go from layer to layer 
(Fig§. 

With a barrier the motion is somewhat more chaotic, 
and the local drift velocity adjusts better to X, as seen in 
Fig[3^. A strict adiabatic approximation would be applica- 
ble if the momentary motion (for a frozen piston position) 
were chaotic, with some finite correlation time t c \. Then 
the adiabatic condition would be Qr c i <C 1. This would 
require to consider a 2D ring, say a Sinai billiard. Within 



this approximation we could define TZ(ip) such that Eq.( 10 ) 
would lead to an Q free parametric integral for Q. 

Contrasting with Boltzmann. — There are two 
conspicuous differences between our results and the ex- 
pectations on the basis of the Boltzmann picture, (a) In 
the Boltzmann picture, with "high" piston, in the ab- 
sence of a barrier, we expect parametric transport with 
1Z = 1, hence obtaining Q = upon integration, im- 
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plying zero drift velocity. (b) Including a barrier, the 
Boltzmann picture suggests a non-zero net transport in 
the same direction that is implied by the pumping oper- 
ation, with < TZ((p) < 1. This means u r > for all the 
layers. Both (a) and (b) are in contradiction to what we 
observe. In particular we point out that the negative u r 
characterizing some of the layers is due to the possibility 
of having a free ballistic contribution with branching ratio 
that favors u < motion. 

One should appreciate the essential difference between 
the Boltzmann and adiabatic pictures: Both would agree 
qualitatively if during the time of g p = 1 the drift velocity 
were zero. Instead it remains constant. One realizes that 
there is an "order of limits" issue: the Boltzmann picture 
assumes that there is always some infinitesimal reflection 
that allows in the adiabatic limit a randomization of the 
velocity. In contrast to that in the adiabatic picture the 
reflection during the g p = 1 stage is strictly zero and hence 
(v) remains constant. 

Summary and discussion. — A phase-space based 
approach for the analysis of stirring in a deterministic 
driven system has been presented. Our oscillating-piston 
model exhibits a layered mixed phase-space structure. 
The determination of the drift velocity requires to go be- 
yond a simple parametric theory: in general neither an 
adiabatic nor a Boltzmann picture applies. The drift ve- 
locity in some layers can even have a sign opposite to the 
current direction that would be expected for a strictly adi- 
abatic pumping ("counter stirring"). These chaotic layers 
appear already for slow driving, whereas a homogeneously 
chaotic phase-space, compatible with a stochastic picture 
of stirring, requires a different limit. It is important to re- 
alize that no simple relation can be established between a 
stirring problem and its corresponding pumping problem 
(that is, the same driven potential in an open configura- 
tion). Different paradigms are involved. 

A few words are in order regarding the quantum case 
[l3] : In the quantum adiabatic limit 1Z can be calculated 
using the Kubo formalism. It has a wide distribution, 
which is the "geometric conductance" analogue of univer- 
sal conductance fluctuations (UCF). One can even observe 
a counter-stirring effect (1Z < 0) which would be impossi- 
ble in the strict classical adiabatic picture. The manifesta- 
tion of dynamical localization requires significantly longer 
time scale of coherence, as explained in Ref. [7] with regard 
to Hamiltonian ratchets. 

Appendix: adiabatic trajectories. — In this ap- 
pendix we clarify what are the equations that describe 
adiabatic trajectories for the pertinent two types of mo- 
tion of the Fermi-Ulam "box model" and its "ring model" 
variation. 

Bounded adiabatic motion: Consider a particle that is 
moving back and forth in a Box of length L, between an 
infinite barrier and an oscillating piston. In the nth colli- 
sion the change in its velocity is —2X(t n ). Summing over 
collisions, approximating by an integral over dt/(2L/^), 



one obtains v(t) = const — (vE/L)X(t), where the abso- 
lute value of the velocity ve = (2E/m) 1 ^ 2 is assumed to 
be approximately constant. This implies that the equa- 
tion of the adiabatic curve in the Poincare section for this 
type of motion is 

Ue 

u((f) = const H cos((p) (25) 

Note that in our ring setting i is 1+ or t_ depending 
whether it is u > or u < trajectory. 

Extended adiabatic motion: There is a similar deriva- 
tion for a particle on a clean ring of length L. There 
is no barrier. The particle changes direction each time 
that it collides with the moving piston. Consequently v(t) 
changes sign each period, and it is convenient to sum the 
increments pairwise. This leads, after an even number of 
collisions, to the obvious result v(t) = const + X(t). This 
implies that the equation of the adiabatic curve in the 
Poincare section for this type of motion is 

u((p) = const + sin((/?) (26) 
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